MIPT/F  ALT-1999-010 


European  Office  of  Aerospace  Research  and  Development 
Project  SPC-98-4036  : 

Contract  F6 1 775-98- WE039 


DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
May  1999  Distribution  Unlimited 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 

Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. _ 

1.  AGENCY  USE  ONLY  (Leave  blank)  I  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 


May  1999 


Final  Report 


4.  TITLE  AND  SUBTITLE 


5.  FUNDING  NUMBERS 


Receptivity  of  Hypersonic  Non-Gradient  Boundary  Layer  to  Wall  Disturbances 


F61775-98-WE039 


6.  AUTHOR(S) 


Dr.  Alexander  Fedorov 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Moscow  Inst,  of  Physics  and  Technology 
16  Gagarin  St, 

Zhukovsky  140160 
Russia 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
EOARD 

PSC  802  BOX  14 
FPO  09499-0200 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

SPC  98-4036 


11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTFtACT  (Maximum  200  words) 


This  report  results  from  a  contract  tasking  Moscow  Inst,  of  Physics  and  Technology  as  follows:  The  contractor  will  perform  theoretical 
analysis  of  the  hypersonic  boundary-layer  receptivity  as  outlined  in  his  proposal  originally  submitted  in  Aug  '97  and  amended  in  Feb  '98. 


14.  SUBJECT  TERMS 


15.  NUMBER  OF  PAGES 


EOARD,  Boundary  Layer  transition,  Hypersonic  Flow 


75 

16.  PRICE  CODE 
N/A 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION  19,  SECURITY  CLASSIFICATION 
OF  THIS  PAGE  OF  ABSTRACT 


20.  LIMITATION  OF  ABSTRACT 


_ UNCLASSIFIED 

NSN  7540-01-280-5500 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  239-18 
298-102 


Abstract 


Theoretical  analysis  of  hypersonic  boundary-layer  receptivity  to  wall  disturbances  is 
conducted  using  a  combination  of  asymptotic  and  numerical  methods.  Excitation  of 
the  second-mode  waves  by  distributed  and  local  forcing  on  the  flat  plate  surface  is 
studied  under  adiabatic  and  cooled  wall  conditions.  Analysis  addresses  receptivity  to 
wall  vibrations,  periodic  suction-blowing  through  a  hole  or  slot,  and  temperature 
disturbances.  A  strong  excitation  occurs  in  local  regions  where  forcing  is  in  resonance 
with  normal  waves.  It  is  revealed  that  the  receptivity  function  tends  to  infinity  as  the 
resonance  point  tends  to  the  branch  point  of  discrete  spectrum  that  is  typical  for  the 
cooled  wall  case.  Asymptotic  analysis  resolves  this  singularity  and  provides  maximal 
receptivity  levels  in  the  branch-point  vicinity.  Analytical  results  are  integrated  into  the 
computational  module  providing  fast  estimates  of  receptivity  levels  for  different  types 
and  shapes  of  wall  forcing.  Numerical  results  indicate  extremely  high  receptivity  to 
vibrations  and  suction-blowing  near  the  lower  neutral  branch.  Critical  amplitudes  of 
local  and  distributed  vibrations  are  estimated  for  bypass  of  the  linear  stability  phase. 
The  theoretical  model  can  be  used  to  predict  initial  amplitudes  of  unstable  waves  in 
active  (ascending)  flights  accompanied  by  skin  vibrations. 
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Nomenclature 


Latin  symbols 

A 

A 

c 

F  =  avJU] 2 

G 

g 

L 

M 

N 

P 

Q 

q 

R  -  4u*ex*  /  v* 

t 

T 

u 

V 

u,v,w 

x,y,z 

JC,  =  X*  /  L  =  £X 


eigenfunction  vector  (direct  problem) 

eigenfunction  vector  (conjugate  problem) 

phase  speed 

frequency  parameter 

receptivity  function 

forcing  shape-function 

streamwise  scale 

Mach  number 

amplification  factor 

pressure  disturbance 

mass-flux  disturbance 

receptivity  factor 

Reynolds  number 
time 

mean  flow  temperature 
mean  flow  streamwise  velocity 
mean  flow  vertical  velocity 
velocity  disturbance 
Cartesian  coordinates 

slow  streamwise  variable 


Greek  symbols 

e-8 1  L 
A 

S  =  Jv\LIU) 


wave-vector  components  in  x-  and  z-direction  respectively 
small  parameter 
delta  function 

boundary  layer  thickness  scale 
kinematic  viscosity 
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pi.cc,  /?) 

Fourier  component  of  the  shape-function  g(x',z') 

CO 

angular  frequency 

e 

temperature  disturbance 

'F 

disturbance  vector-function 

F(x,y) 

disturbance  amplitude  vector-function 

resonance  point 
adiabatic  wall 
branch  point 

upper  boundary-layer  edge 

imaginary  part 

real  part 

suction-blowing 

vibrations 

wall 

temperature  perturbation 


Superscripts 

bar 

complex  conjugate  value 

asterisk 

dimensional  value 

Abbreviations 

PSE 

parabolized  stability  equations 

Im 

imaginary  part 

Re 

real  part 
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1.  Introduction 


In  hypersonic  flights  and  quiet  wind  tunnels,  freestream  perturbations  and  body- 
induced  disturbances  are  normally  small.  In  this  case,  the  transition  process  can  be 
treated  as  an  initial  boundary  value  problem,  which  comprises  excitation  of  unstable 
boundary-layer  modes  (receptivity  problem),  their  downstream  amplification  (linear 
stability  problem)  and  nonlinear  breakdown  to  turbulence  [1,2].  Transition  prediction 
technology  is  usually  based  on  the  eN  -method  [3-5],  which  considers  only  the  linear 
phase  of  instability  growth.  All  other  physics  including  receptivity  and  nonlinear 
breakdown  is  incorporated  into  the  value  N,  which  is  colibrated  for  a  particular  set  of 
experimental  data.  The  finite  amplitude  method  including  receptivity,  instability 
amplification,  and  an  amplitude  criterion  for  transition  onset  can  avoid  many 
shortcomings  of  the  eN  -method.  Recent  codes  based  on  the  parabolized  stability 
equations  (PSE)  can  predict  the  instability  evolution  and  capture  the  initial  phase  of 
nonlinear  breakdown  [6-8].  However  the  PSE  codes  require  initial  distributions  of  the 
disturbance  amplitude  as  an  input,  which  should  be  supplied  by  receptivity  theory. 

Despite  of  a  good  progress  in  experimental  studies  of  hypersonic  stability  [9],  we 
found  only  few  papers  relevant  to  experiments  on  hypersonic  receptivity  [10,  11]. 
Kendall  [10]  measured  correlation  between  the  freestream  and  boundary-layer 
disturbances  in  the  Mach  number  range  1. 6-5.6,  and  found  that  the  correlation 
coefficient  is  growing  with  the  Mach  number.  Experimental  data  obtained  on  a  flat 
plate  at  Mach  4.5  are  in  a  good  agreement  with  the  forcing  theory  of  Mack  [12]. 
However  Mack  noted  that  the  major  difficulty  in  the  use  of  the  forcing  theory  is  that 
force  disturbances  are  distinct  from  free  disturbances,  and  the  process  by  which  the 
former  becomes  the  latter  is  unknown.  Maslov,  Shiplyuk,  Sidorenko  and  Amal  [11] 
studied  the  flat-plate  leading  edge  receptivity  to  acoustic  disturbances  in  the  Mach  6 
wind  tunnel.  Acoustic  field  was  generated  by  a  point  source  and  plane  source.  It  was 
found  that  the  leading  edge  receptivity  depends  on  the  acoustic  wave  inclination  angle 
X .  Maximum  receptivity  was  observed  at  x  *  60° .  The  receptivity  coefficient  for 
two-dimensional  acoustic  waves  with  the  phase  speed  e  =  1  - 1  /  Me  was  much  larger 
than  that  measured  for  three-dimensional  acoustic  waves  generated  by  the  point 


10 


source.  These  observations  are  consistent  with  the  theoretical  results  of  Fedorov  and 
Khokhlov  [13]. 

Only  few  theoretical  studies  have  been  addressed  hypersonic  receptivity  problem. 
Fedorov  and  Khokhlov  [13]  analyzed  receptivity  on  the  sharp  leading  edge  of  a  flat 
plate  to  external  disturbances  including  sound  waves,  entropy  and  vortical  waves  as 
well  as  leading-edge  vibrations.  They  showed  that  freestream  noise  and/or  vortical 
disturbances  induce  a  local  acoustic  source  on  the  leading  edge.  This  source  generates 
acoustic  waves,  which  are  synchronized  with  boundary-layer  normal  modes  (the  first 
and  second  modes  according  to  Mack’s  classification).  That  causes  a  strong  excitation 
of  unstable  normal  waves  by  freestream  disturbances  through  the  mechanism,  which 
is  qualitatively  different  from  that  observed  at  subsonic  speeds  [14-19].  Zhong  [20] 
performed  direct  numerical  simulation  of  the  parabolic  leading-edge  receptivity  to 
freestream  sound  in  a  hypersonic  flow.  In  his  calculations,  shock-layer  modes  and 
unstable  boundary-layer  modes  were  detected  between  a  bow  shock  and  body  surface. 
Fedorov  [21]  showed  that  external  vorticity  and  entropy  disturbances  can  effectively 
generate  the  second  mode  near  the  lower  neutral  branch.  In  this  region,  the  first  and 
second  mode  waves  can  be  synchronized  with  entropy/vorticity  waves  of  the  phase 
speed  c  =  1 .  In  the  synchronism  region,  entropy/vorticity  disturbances  are  partially 
swallowed  by  the  boundary  layer  and  excite  instability.  This  mechanism  can  compete 
with  the  leading-edge  receptivity  in  cases  of  conical  configurations  such  as  a  sharp 
cone. 

Another  important  type  of  disturbances  is  associated  with  body-induced  perturbations. 
In  active  flights,  propulsion  system  generates  vibrations  of  the  vehicle  skin.  These 
vibrations  can  propagate  upstream  and  excite  the  boundary-layer  instability  on  the 
forebody  surface.  Skuratov,  Fedorov  and  Shogin  [22]  published  transition  data 
obtained  on  a  7.5°  half-angle  sharp  cone  installed  on  the  rocket  nose.  Transition  loci 
were  measured  during  ascending  flights  in  the  local  Mach  number  range  2.5-4  at 
wall  temperature  ratios  TwITad  =  0.2 -0.8.  Transition  onset  points  were  close  to 
those  obtained  in  conventional  (noisy)  wind  tunnels.  Such  an  early  transition  seemed 
to  be  due  to  vibrations  of  the  cone  surface.  Another  source  of  disturbances  is 
associated  with  gaps  and  slots  between  elements  of  the  thermal  protection  system. 
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They  may  work  as  resonators,  which  induce  periodic  suction-blowing  on  the  body 
surface.  These  perturbations  can  also  generate  unstable  modes  in  the  boundary  layer. 

Theoretical  and  experimental  studies  of  the  boundary-layer  receptivity  to  body- 
induced  perturbations  have  been  performed  for  subsonic  flows.  It  was  found  that  the 
boundary  layer  is  extremely  sensitive  to  wall  perturbations  to  be  in  resonance  with 
Tollmien-Schlichting  waves.  For  example,  the  theoretical  analysis  [23]  and  the  wind- 
tunnel  experiment  [24]  showed  that  two-dimensional  vibrations  of  5  microns 
amplitude  (about  1%  of  the  boundary-layer  thickness)  generate  the  TS  waves  with  the 
initial  amplitude  Ain  &  0.1%.  Assuming  that  nonlinear  breakdown  occurs  at  the 
critical  amplitude  Acr  « 1% ,  we  obtain  that  the  amplification  factor  is  N  = 
In (Acr  IAin )  «2.3 .  This  value  is  much  less  than  the  empirical  value  N  «  10  typical 

for  “quiet”  conditions.  Similar  conclusion  has  been  made  regarding  periodic  suction¬ 
blowing  induced  through  a  two-dimensional  slot  in  subsonic  flat  plate  [25].  These 
findings  lead  to  the  assumption  that  hypersonic  boundary  layers  may  be  also  very 
receptive  to  wall  perturbations. 

In  this  report,  we  perform  theoretical  analysis  of  hypersonic  boundary-layer 
receptivity  to  wall  disturbances.  We  consider  excitation  of  the  second-mode  waves  by 
distributed  and  local  forcing  on  a  flat  plate.  In  Section  2  we  formulate  the  problem, 
analyze  disturbance  field  in  the  case  of  resonance  between  the  second-mode  wave  and 
distributed  wall  forcing.  Then  we  consider  receptivity  to  local  forcing  produced  by  a 
vibrator,  periodic  suction-blowing  through  a  slot  or  hole,  and  temperature 
perturbations  induced  by  a  heating  element.  In  Section  3  we  describe  computational 
algorithm  providing  the  receptivity  function  and  discuss  numerical  results  for 
adiabatic  and  cooled  plate  cases.  We  show  that  the  receptivity  function  tends  to 
infinity  as  the  resonance  point  tends  to  the  branch  point  of  the  discrete  spectrum. 
Receptivity  to  wall  forcing  in  the  branch-point  vicinity  is  analyzed  in  Section  4. 
Critical  amplitudes  of  wall  vibrations  providing  bypass  of  the  linear  stability  phase  are 
estimated  in  Sections  5  for  local  and  distributed  forcing.  In  Section  6,  we  conclude  the 
report  and  discuss  future  effort. 
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2.  Receptivity  to  wall  perturbations  in  the  case  of  a  simple 
discrete  spectrum 

In  this  section,  we  analyze  the  disturbance  field  generated  by  wall  perturbations  in  the 
case  of  a  simple  discrete  spectrum;  i.e.  eigenvalues  of  the  boundary-layer  modes  are 
assumed  to  be  different.  Under  this  condition,  analysis  of  hypersonic  boundary-layer 
receptivity  to  wall  forcing  can  be  performed  using  the  theoretical  model  of  [23,  25,  26] 
developed  for  subsonic  boundary  layers. 

2.1  Problem  formulation  for  distributed  forcing 

We  consider  a  two-dimensional  laminar  boundary  layer  of  perfect  gas  on  a  hypersonic  flat 
plate  schematically  shown  in  Fig.  1 .  It  is  assumed  that  the  global  streamwise  length  L  is 

much  larger  than  the  boundary  layer  thickness  scale  8  =  -Jv’L  / U]  ;  i.e.  the  ratio 
s  —  8  /  L  is  small.  This  is  equivalent  to  the  assumption  that  the  Reynolds  number 
RL=Juiu7,=e-'  is  large.  The  longitudinal  x* ,  normal  y* ,  transversal  z* 
coordinates  and  time  t*  are  made  nondimensional  using  8  and  U*e  as 

(x,y,z)  =  (x\y*,z*)/S,  t  =  t*U*e  IS. 

We  introduce  the  slow  variable  x,  =  x*  /  L  =  ex  and  specify  the  nondimensional  mean- 
flow  velocity  components  (U,  V)  and  temperature  T  as 


u^u'/u:=  t/(x,  ,y),  V  =  V*  /  u:  =  £V0(XJ  ,y),  T  =T*  /  TJ  =  T(x,  ,y) .  (2.1) 

Distributed  perturbations  of  velocity  (u,v,w),  pressure  p  and  temperature  0  are  induced 
by  vibrations  or  any  other  forcing  on  the  body  surface.  Wall  perturbations  are  represented 
in  the  traveling- wave  form 

(u,v,w,0)  =  (q>ltg>2,q>3iq>4)exp(iooc  +  ifiz-ia)t)  at  y  =  0,  (2.2) 
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where  a  and  /?  are  wavenumber  components;  <±>  is  angular  frequency.  The  forcing 
amplitude  is  assumed  to  be  small  and  the  following  analysis  is  performed  using  a  linear 
theory  of  small  perturbations.  The  flow  disturbance  is  described  by  the  vector  function 


du  n  39 

vF  =  (w,— ,v,p,0,— ,w,— )  , 
dy  dy  dy 


T( x,y,z,t )  =  F(x,y)exp(ifiz-iat) . 


(2.3) 

(2.4) 


The  amplitude  function  F(x,y)  satisfies  a  system  of  partial  differential  equations 
resulted  from  Fourier  transform  of  linearized  Navier-Stokes  equations  with  respect  to 
time  and  z  -coordinate.  These  equations  can  be  written  in  the  matrix-operator  form 

H(y,  dy ,  x, ,  sdXi ,  co,  fi)F  =  0 .  (2.5) 


Boundary  conditions  are  specified  as 


(Fx,F3,Fs,F1)  =  {(px,cp2,(pi,(pi) exp(/ ax)  at_y  =  0,  (2.6) 

|f|  <  oo  aty  — >  oo .  (2.7) 

As  contrasted  to  a  standard  stability  analysis,  the  problem  (2.5)-(2.7)  has  non- 
homogeneous  boundary  conditions  on  the  wall  (y  =  0)  and  describes  the  flow  disturbance 
field,  which,  in  general,  includes  forcing  and  normal  waves  generated  by  wall 
perturbations. 

2.2  Non-resonance  case 

We  consider  a  partial  solution  of  the  problem  (2.5)-(2.7)  as  a  decomposition  of  the 
discrete  modes  Fk(xx,y)  with  eigenvalues  ak{xx,p,co)  and  the  force  term  Fv(x,,y) 
with  the  wavenumber  a .  Then  the  disturbance  amplitude  is  expressed  in  the  form 
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(2.8) 


^  =  £(4,  +  sFkl+---)exp(ie  lSk)  +  Fv  exp  (is  'ax,), 

k 

sk  =  \ak(xl)dxl.  (2.9) 


Substituting  (2.8)  into  Eqs.  (2.5)-(2.7)  and  grouping  terms  of  the  same  order  of  magnitude 
with  respect  to  s  we  obtain  a  sequence  of  problems  for  the  amplitude  functions 
Fkj  (x,  ,y) ,  j  =  0,  1 ,  ...  In  the  first  order  approximation,  we  get  the  eigenvalue  problem 


f  d_ 
\dy 


-Hn 


FkO  ~  0  5 


(-^lo)  1.3,5. 7  —  0  aty  —  0, 

\fJ  — >  0  at y  ->  oo. 


(2.10a) 

(2.10b) 

(2.10c) 


Here  the  matrix  H0  has  dimension  8  x  8 ;  its  elements  depend  on  the  mean-flow  profiles 
U{xx,y)  and  T(xx,y),  disturbance  parameters  a,  (3,  to  and  Reynolds  number 

R  =  *ju*x*  /  v*e  .  Explicit  form  of  the  matrix  H0  is  given  in  [26,  27].  Equations  (2.10a)- 
(2.10c)  form  a  standard  linear  stability  equation  system  for  locally  parallel  mean  flows. 

Forcing  term  Fy  (x,  ,y)  is  a  solution  of  the  non-homogeneous  problem 


f 

v 


^  ° 


=  0, 


(K  ) I.3.5.7  —  ^1,2,3.4  y  ~  0  , 


(2.11a) 

(2.11b) 


<  oo  at  y  -»  oo . 


(2.11c) 


If  the  forcing  wavenumber  (a,  (3)  coincides  with  a  wavenumber  of  continuous  spectrum, 
then  the  body  surface  effectively  generates  acoustic,  vortical  or  entropy  waves.  If  {a,  ft) 
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does  not  belong  to  the  continuous  spectrum,  then  the  forcing  disturbance  decays 
exponentially  as  y  -»  oo .  Hereafter  we  analyze  the  latter  case. 

In  the  second  order  approximation,  we  obtain  the  non-homogeneous  problem 


dHn  dFm 

(2.12a) 

H,)FkX 

=  "  da  A,  +H'Fu ’ 

iFkX\ 

,3,5,7  =  0  aty  =  0 , 

(2.12b) 

|4| 

->  0  at  y  -»  oo . 

(2.12c) 

An  explicit  form  of  the  matrix  //,  is  given  in  [26,  27],  Its  elements  depend  on  the  mean- 
flow  profile  derivatives  dU  /  dct,  oT  /  dcx  and  vertical  velocity  V0.  This  matrix  is 
associated  with  nonparallel  effects  due  to  the  boundary-layer  downstream  growth. 

Solution  of  the  eigenvalue  problem  (2.10)  can  be  represented  in  the  form 

Fk0  =  ck(xx')Ak(xx,y,a k) ,  (2.13) 

where  Ak  is  eigenfunction  normalized  by  a  certain  condition.  For  example,  the  pressure 
disturbance  amplitude  is  constant  on  the  wall:  Aki(x{, 0,ak)  =  1 .  Problem  (2.12)  has  a 

nontrivial  solution  if  its  right-hand  side  is  orthogonal  to  the  eigenfunction  Bk{xx,y,ak) 
of  the  conjugate  problem.  This  leads  to  the  following  ordinary  differential  equation  for 
the  amplitude  coefficient  ck  (x, ) 


-  dH0  -  dck  r  •  n 

i  <  Bk,—±Ak  >  —  =  [~i  <  Bk 

da  dxl  da  dc 


=  m°dAt  >+<Bl,HtAt>¥^ 


ck(xu)  =  ch. 


(2.14) 

(2.15) 
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where  the  scalar  product  is  defined  as 


CO  g 


<B,A>=feBJAJdy 


(2.16) 


0  >1 


Upper  bar  denotes  complex  conjugate  values;  the  amplitude  coefficient  ck  is  assumed  to 
be  prescribed  at  some  initial  point  . 


If  the  interaction  between  normal  modes  is  neglected,  then  Eq.  (2.14)  describes  the 
normal  wave  propagation  in  a  weakly  nonparallel  boundary  layer.  This  approximation  is 
valid  in  a  region  where  eigenvalues  are  different;  i.e.  ak(xx)  *  a  j(xx)  at  k  *  j ,  and  the 

forcing  wavenumber  a*  ak;  i.e.  in  the  absence  of  resonance  between  wall  forcing  and 
normal  wave,  receptivity  is  weak. 


2.3  Resonance  between  forcing  and  normal  waves 

We  assume  that  the  force  wave  is  in  resonance  with  a  normal  wave  of  discrete  spectrum 
at  some  point  xx  =  x]0;  i.e.  a  »  an(xx0)  =  a0 .  We  consider  the  case  when  the  normal 
wave  is  unstable  and  the  resonance  point  x]0  is  close  to  the  lower  neutral  branch.  It  is 
expected  that  the  wall  disturbance  generates  the  normal  wave  in  a  local  region  near  the 
resonance  point.  Downstream  from  this  region,  the  normal-wave  amplitude  exponentially 
increases  due  to  the  boundary  layer  instability. 

Flow  scheme  of  resonant  forcing  is  shown  in  Fig.  2.  Analysis  of  Eqs.  (2.5)-(2.7)  indicates 
that  a  strong  excitation  occurs  in  a  small  region  where  |x,  -  x]0|  =  0(s112 ) .  We  introduce 
the  inner  longitudinal  variable  as 

Z=£-]/\xx-xx0)  =  O(l).  (2.17) 
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In  dimensional  form,  the  resonance  point  is  x*0=L,  and  the  inner-region  length  is  of  the 
order  of  /  =  L  •  Rj'12  =  8  ■  R^2 .  This  length  is  much  larger  than  the  boundary-layer 
thickness  scale  8  and  much  smaller  than  the  distance  L  from  the  plate  leading  edge  to 
the  resonance  point. 

We  assume  that  the  force  wavenumber  is  a  =  aQ+ysV2 ,  where  the  real  parameter  y 
characterizes  resonance  detuning.  Then  the  inner  solution  of  Eqs.  (2.5)-(2.7)  is  expressed 
in  the  form 


F  =  [s~U2c(^)A0(y)  +  Av(£,y)  +  ...]exp  {is~ula^),  (2.18) 


where  A0(y)  =  An(x]0,y,a0)  is  the  n th-mode  eigenfunction  at  the  resonance  point. 
Expanding  the  operator  H0  in  the  vicinity  of  x  =  xl0  we  obtain  the  following  inner 

problem  for  the  force  wave  amplitude  Ar 


H0(x]0,a0) 


A,  =  -i 


dc  dH0 
d%  da 


(xl0,a0)A0  +  c£  (xl0,aQ)A0, 
ox. 


(Ay) 

1,3,5, 7  =  ^.,2,3,4  (^10  )exp(z>£)  at  y  =  0 , 


-»  0  at  y  -»■  oo . 


(2.19a) 

(2.19b) 

(2.19c) 


Using  a  standard  approach  (see  for  example  [28])  we  can  express  the  problem  (2.19)  in 
the  form 


Ho(xl0,ao) 


A.,  =  -i 


^^-(r10,a0H  +  ^U0,ffoH  +  A(#o  exp (iyg), 
dq  doc  (dx,  j 


(2.20a) 

(^v)i,3,5,7  =  °  aty  =  0 , 

(2.20b) 

\x  ->  Oaty->  oo  , 

(2.20c) 
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^0  =  iSPx  (*io  \^(p2  (^10 )» 0,  (p 3  (x10 ),  0,  ^)4  (x10 ),  0) T , 


(2.20d) 


where  A(y)  is  delta  function.  Problem  (2.20)  has  a  non-trivial  solution,  if  its  right-hand 
side  is  orthogonal  to  the  eigenfunction  B0  =  Bn(xX0,y,a0 )  of  the  conjugate  problem. 
This  condition  leads  to  the  following  equation  for  the  amplitude  coefficient  c(E,) 


—  -ib%c  =  qmp(iy%). 

(2.21a) 

b~  "(*io)> 

dxx 

(2.21b) 

(B0,®o)y=0 

q  ~ 1  -  cH 

<-  B0,  (xxo,a0)A0  > 

(2.21c) 

If  the  normal  wave  has  zero  amplitude  upstream  from  the  resonance 

boundary  condition  for  Eq.  (2.21a)  is  specified  as 

region,  then  the 

c(g)  — >  0  at  B,  — >  -oo . 

(2.22) 

Solution  of  the  problem  (2.21a),  (2.22)  is  expressed  in  the  form 

,  f  ibs2) 

c(£)  =  q  exp (ib%  /  2)  1  exp  iys  -  ds . 

J  v  2 

-OO 

(2.23) 

In  typical  cases,  the  unstable  mode  eigenvalue  satisfies  the  following  conditions  near  the 
lower  neutral  branch: 


dRe(an)  /  dxx  >0,  dlm(an)l  dxx  <0. 
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Then  an  asymptotic  behavior  of  the  amplitude  coefficient  c(#)  as  #  +°°  is  associated 

with  the  saddle  point  s,=y/b  and  evaluated  using  the  steepest  descent  method  as 


2m 

C(£)  ->  exp 


^  iy 2  ib?\ 

{2b+  2  , 


at#- 


+CO . 


(2.24) 


Analyzing  asymptotic  trends  of  the  outer  solution  Fn0  we  conclude  that  its  inner  limit  at 
x,  -»  jc10  +0  matches  with  the  outer  limit  of  the  inner  solution  at  #-»  +co ;  z'.c.  in  the 
resonance  region,  x,  -  x10  =  0(evl ) ,  the  force  wave  excites  the  unstable  normal  wave  of 
the  amplitude 


Fn=s'mq 


2tu 


iy 2% 

4  “MF J 


A„(xl,y)exp 


is"'  \an(xx)dxx 


*10 


for  x,  >  x10 .  (2.25) 


This  equation  shows  that  the  normal  wave  amplitude  is  proportional  to  ( L  /  S)xn  =  Rx[2 , 
whereas  the  resonant  region  length  is  /  =  R]/2S .  As  the  Reynolds  number  increases, 
nonparallel  effects  (which  disturb  the  resonance  condition)  decrease  and  the  relative 
length  /  /  8  increases.  In  the  parallel  mean-flow  limit,  RL->  °o ,  the  resonance  region  is 
infinitely  large  and  the  normal  wave  amplitude  tends  to  infinity.  In  this  case,  the  upper 
limit  of  disturbance  amplitude  is  determined  by  nonlinear  effects.  Both  nonlinear  and 
nonparallel  mechanisms  of  the  resonance  detuning  can  be  observed  simultaneously  at 
sufficiently  large  Reynolds  numbers  and  strong  forcing.  Note  that  the  expression  (2.25)  is 
similar  to  that  obtained  by  Tumin  and  Fedorov  [23]  for  subsonic  boundary  layers  using 
another  approach. 


2.4  Local  forcing 

We  consider  a  force  disturbance,  which  is  localized  in  a  region  comparable  with  the 
boundary  layer  thickness.  This  may  be  a  vibrator,  suction/blowing  induced  through  a  hole 
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(in  three-dimensional  case)  or  slot  (in  two-dimensional  case),  and  temperature 
perturbations  induced  by  a  heating  element.  In  the  general  case,  the  force  shape  is 
expressed  in  the  vector  form 

V  =  (u,v,e,wf\y=Q  =  g(x',z')p0,  x'  =  x-x0  =  £~\x{ -X]0),  z'  =  z-z0,  (2.26) 

where  %  is  force  amplitude  vector,  which  does  not  depend  on  the  local  coordinates  x' 
and  z' .  However  this  vector  is  a  function  of  the  force  type  and  mean- flow  characteristics 
at  the  central  point  (x0,z0).  Its  explicit  form  will  be  given  in  Section  3. 

The  force  shape  vector-function  is  associated  with  the  Fourier  transform 


-r<-u 

p(cc,{3)  =  Jdz'  Jdx'g(x',z')exp(-iax'  -i/k'). 


(2.27) 


and  can  be  expressed  as 


_fo_ 

(2n) 


<P  = 


+oo  +oo 

,/?)exp(/ax'  +  i/3z')da  = 


+oo  +oo 


T£U2eia°x'  jdJ3  jp(a0  +  £U2y,/3)exp(iy%  +  iflz')dy . 

-CO  -00 


(2.28) 


For  two-dimensional  forcing  of  the  shape  g(x'),  Equations  (2.27)  and  (2.28)  are  reduced 
to  the  form 


tw 

p(a)=  ^dx'g{x')txp(-iax'). 


(2.29) 


—  +oo 

<P  =  P~  \p{a)txp(iax')da  =  ^-£{l2eia,,s'  \p{a0  +  £my)exp(iy%)dy .  (2.30) 

/  TT  J  /  77"  J 


%  -1/2 


In 


2  n 
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Using  Eq.  (2.25)  we  evaluate  the  normal-wave  amplitude  generated  by  force  waves  in  the 
wavenumber  range  (a0  +  <?'/2(y  +  dy),p  +  dp)  as 


dF„  =■ 


{l7C)2 


2m 

' iy2' 

hexp 

[lb) 

(  A' 


'\p(a0  +£  y,P)An(x',y;p)ex  p 


i 

\  0 


J an(x')dx ' 


eip:  dydfi . 


(2.31) 


Integrating  over  y  we  get  the  following  wave  packet  in  the  first  order  approximation 


2  n 


hoo  a 

\qp(a, ,P)An ( x',y-p ) exp  i \an{x')dx'  +  ipz' 


\  0 


dp. 


(2.32) 


For  two-dimensional  forcing  specified  by  Eqs.  (2.29),  (2.30),  the  normal  wave  amplitude 
is  expressed  as 


(2.33) 


According  to  Eq.  (2.21c)  the  receptivity  factor  q(xw,a0 )  does  not  depend  on  the  force 
shape  g(x',z').  If  this  factor  is  calculated  once  and  for  all,  then  we  can  evaluate  the 
normal  wave  amplitude  induced  by  local  forcing  of  any  shape. 

As  contrasted  to  the  case  of  distributed  forcing  (see  Eq.  (2.25)),  local  forcing  generates 
the  boundary-layer  normal  waves  of  a  finite  amplitude  even  in  the  parallel  mean  flow 
limit.  This  is  due  to  the  fact  that  only  a  small  portion  of  the  force  spectrum, 
(a  -  a0)  =  0(sm),  is  involved  into  the  resonance  mechanism.  As  s  -»  0,  the  resonant 
excitation  increases  proportionally  to  s~m ,  whereas  the  spectrum  region  is  narrowed 
down  proportionally  to  sm .  Because  both  trends  compensate  each  other,  the  normal- 
wave  amplitude  excited  by  local  forcing  tends  to  a  finite  value. 
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3.  Numerical  results 


3.1  Computational  algorithm  for  receptivity  to  local  forcing 

Summarizing  results  of  Section  2  we  formulate  the  following  algorithm  for  calculations 
of  the  normal-wave  amplitude  generated  by  local  wall  forcing: 

1.  Specify  the  function  <p(x',z')  defined  by  Eq.  (2.26),  where  the  vector  <jz>0  depends  on 
a  forcing  type.  For  vibrations,  it  is  expressed  as 

V*  =  (U\  ieo,  r,  0)^O  V]=V|o ,  (3.1) 

where  prime  denotes  the  partial  derifative  d  I  dy .  In  this  case,  the  force  shape  function 
g(x',z')  corresponds  to  the  wall  displacement  amplitude. 

For  suction-blowing  inducing  a  vertical  velocity  on  the  body  surface,  we  obtain 

<p0=(0,l,0,0)r.  (3.2) 

In  this  case,  the  function  g(x',zr)  determines  the  vertical  velocity  amplitude. 

For  temperature  perturbations  on  the  wall  surface,  we  get 

^0=(0,0,l,0)r.  (3.3) 

In  this  case,  the  function  g(x',z')  determines  the  surface  temperature  amplitude. 

2.  Solve  direct  and  conjugate  eigenvalue  problems  at  the  central  point  xl  =  xl0  of  the 
force  element  and  calculate  the  normal- wave  eigenvalue  aQ  =  an(xlQ,fi,oj)  and 
eigenfunctions  A0(y) ,  B0(y) .  This  can  be  done  using  a  standard  stability  code. 
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3.  Calculate  the  receptivity  factor  q(x]0,a0,J3,co)  using  Eq.  (2.21c). 

4.  Evaluate  the  normal-wave  amplitude  Fn  using  Eq.  (2.32)  for  three-dimensional 
forcing  or  Eq.  (2.33)  for  two-dimensional  forcing. 

This  algorithm  was  integrated  in  a  computationally  non-intensive  module,  which  was 
coupled  with  a  standard  stability  code.  This  module  works  as  a  black  box  and  predicts 
initial  amplitudes  of  unstable  normal  waves  induced  by  wall  forcing  of  prescribed  shape. 

3.2  Numerical  results 

Calculations  have  been  conducted  for  the  Mack  second-mode  induced  by  a  two- 
dimensional  local  forcing.  For  surface  vibrations,  the  shape  vector-function  <pb  is 
determined  from  Eq.  (3.1).  The  disturbance  amplitude  is  characterized  by  the  mass-flux 
perturbation  evaluated  at  the  point  of  its  maximum:  y  -  ym .  This  quantity  is  expressed  as 

Q(x]0,ym)=\  p(a0)\Gvih,  (3.4) 

where  Fourier  component  p{a)  is  determined  from  Eq.  (2.29);  the  receptivity  function 
for  local  vibrations  is  calculated  using  Eqs.  (2.20d),  (2.21c),  (2.33)  and  (3.1)  as 


Gvjb(X  10 ’^0  ) : 


02  O'*.  *10  > 


(U  B2  ]  +  ICoB 2  3  +  T B2  5  )  v=0„v,  =.v„ 


i  da 


(•*•10  >a0  )-^2  ’^2 


(3.5) 


This  function  characterizes  the  mass-flux  amplitude  of  the  normal  wave  generated  by  a 
vibrator  with  the  unit  Fourier  component  p(a0 )  =  1 ;  Q2  (y„, ,  xl0 )  is  mass  flux  disturbance 

evaluated  using  the  second-mode  eigenfunction  A2.  For  the  boundary  layer  on  a  flat 

plate,  the  length  scale  is  determined  as  8  =  ^ v'ex’0  /  U'e  ,  and  the  Reynolds  number  is 
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R  =  ^U*ex o  /  ve  .  If  the  receptivity  function  Gvjh  is  known,  then  we  can  calculate  the 

initial  amplitude  of  the  second-mode  wave  generated  by  a  vibrator  of  arbitrary  shape 
using  Eqs.  (3.4)  and  (2.33). 

In  a  similar  way,  we  evaluate  receptivity  to  suction-blowing  through  a  two-dimensional 
slot.  The  mass-flux  perturbation  is  expressed  in  the  form 


Q(xw,ym)=\p{a^)\Gs,  (3.6) 

where  the  receptivity  function  is  determined  from  Eqs.  (2.20d),  (2.21c),  (2.33)  and  (3.2) 
as 


Gs(x  ,o?«o)  = 


Q2(y,n^0)- 


)v=0,.v,-\„ 


dju 

da 


(*10  ,«o  )^2  ’^2 


(3.7) 


The  mass-flux  disturbance  excited  by  surface  temperature  perturbations  is  expressed  in 
the  form 


Q(X\«,ym)=\p{a,)\Ge’  (3-8) 

where  the  receptivity  function  is  determined  from  Eqs.  (2.20d),  (2.21c),  (2.33)  and  (3.3) 
as 


(x,0,a0 )  ~ 


&  (T,„, *,<>)- 


(■®2.5  )  v=0,.v,=.vlr 


da 


(*10>a0  )^2  >-®2 


(3.9) 
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Case  1:  Cooled  wall 

Calculations  of  the  receptivity  functions  Gvib ,  Gs  and  Ge  versus  the  Reynolds  number 
R  -  -Ju'exl  /  v*  have  been  made  for  the  following  mean  flow  parameters:  local  Mach 
number  Me  =  5.95 ,  ffeestream  temperature  T*  =  70.2578  K,  wall  temperature  ratio 
r;  /  T'e  =  0.69539  corresponding  to  the  ratio  Tw  /  Tad  =  0.1 ,  specific  heat  ratio 
Cp  /  cv  =  1.4 ,  Prandtl  number  Pr=0.72,  second  viscosity  coefficient  k  =  1.2 .  Disturbance 

frequency  parameter  is  F  =  10'4.  As  shown  in  [21]  the  discrete  spectrum  has  a  singular 
topology  in  this  case.  First  and  second  modes  have  the  branch  points  Rb]  and  Rb2  located 
near  the  lower  and  upper  neutral  branch  respectively.  The  first  branch  point  is  very  close 
to  the  real  axis  of  complex  R  -plane.  The  second  branch  point  is  slightly  shifted  from  the 
real  axis.  In  the  vicinity  of  the  branch  points,  the  eigenvalue  spectrum  splits  as  shown  in 
Figs.  3a  and  3b  (see  also  Figs.  14a  of  Ref.  [21]  for  a  conical  boundary  layer). 

Figures  4a,  4b,  4c  and  4d  show  the  second-mode  wave  increment  -a,  =  -lma2(R)  and 
the  receptivity  function  Gvih{R),  GS(R)  and  Ge(R)  respectively.  It  is  seen  that  the 
receptivity  functions  have  a  sharp  peak  at  the  Reynolds  number  R  «  Rbi.  As  shown  in 
[21],  the  scalar  product  has  the  following  trend  near  the  branch  point 

as  =  1,2.  (3.10) 

According  to  Eqs.  (3.5),  (3.7)  and  (3.9)  the  receptivity  functions  tend  to  infinity  as 
R  -»•  Rhl.  Second  maximum  of  Gvjh(R) ,  GS(R)  and  G0(R)  is  observed  at  the  Reynolds 

number  R  «  2500 ,  which  is  close  to  the  second  branch  point  Rbl .  Because  the  point  Rhl 
is  shifted  to  the  complex  R  -plane  this  maximum  is  smaller  and  smoothed  out. 
Comparing  data  in  Figs.  4b-4d  we  conclude  that  receptivity  to  suction-blowing  is  several 
times  higher  than  receptivity  to  vibrations.  In  turn,  receptivity  to  vibrations  is  several 
times  higher  than  receptivity  to  temperature  perturbations. 
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Figures  5a,  5b,  5c  and  5d  show  distributions  of  -a,(R) ,  Gvjh(R),  GS(R )  and  Ge(R)  at 
various  frequency  parameters  F.  The  highest  peak  of  the  receptivity  function  corresponds 
to  Curve  3  at  F  =  10-4 .  In  this  case,  the  first  branch  point  Rbl  is  very  close  to  the  real  axis 

of  the  complex  R  -plane  as  shown  in  Fig.  6.  For  the  frequency  parameter  F  <  10“4 ,  the 
first  branch  point  is  located  in  the  upper  half  of  the  complex  R -plane  (Points  4  and  5  in 
Fig.  6).  As  the  Reynolds  number  increases  along  the  real  axis,  this  branch  point  is 
bypassed  from  the  lower  side.  The  first  mode  (with  the  phase  speed  behavior 
c  — » 1  + 1  /  Me  as  R  decreases)  becomes  unstable.  In  this  case,  the  phase  speed  crosses  the 
ray  c  =  1  that  causes  discontinuity  of  the  receptivity  function.  For  the  frequency 
parameter  F  >  10'4 ,  the  branch  point  is  in  the  lower  half  of  the  complex  i?-plane  (Points 
3,  2  and  1  in  Fig.  6).  As  the  Reynolds  number  increases  along  the  real  axis,  this  branch 
point  is  bypassed  from  the  upper  side.  The  second  mode  (with  the  phase  speed  behavior 
c  — >  1  - 1  /  Me  as  R  decreases)  becomes  unstable.  Note  that  both  cases  correspond  to  the 
second-mode  instability  in  terms  of  Mack’s  classification. 

Figure  7a  shows  that  the  receptivity-function  maximum,  Gribm  -  max[G,.,.6(R)],  is  close 
to  the  lower  neutral  branch.  Figures  7b  and  7c  illustrate  that  the  function  Gvjbm(F )  has  a 

sharp  peak  at  the  frequency  F  =  10"4  corresponding  to  a  singular  case  when  the  branch 
point  crosses  the  real  axis  (see  Point  3  in  Fig.  7c).  The  receptivity  functions  for  suction¬ 
blowing  perturbations,  Gs ,  and  temperature  perturbations,  Ge ,  have  similar  behavior. 

Summarizing  we  conclude  that  hypersonic  boundary-layer  on  a  cooled  plate  is  very 
sensitive  to  the  vertical  velocity  perturbations,  if  they  are  in  resonance  with  boundary- 
layer  modes.  Such  perturbations  may  be  induced  by  suction-blowing  through  a  hole  or  by 
skin  vibrations.  Receptivity  to  temperature  perturbations  is  relatively  low.  The  receptivity 
level  strongly  depends  on  the  forcing  locus.  The  receptivity  functions  have  a  sharp  peak 
in  the  vicinity  of  the  branch  point  Rbl .  This  peak  is  observed  near  the  lower  neutral 
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branch  of  the  second  mode.  Singular  behavior  of  the  receptivity  functions  near  the  point 
Rh]  indicates  that  the  foregoing  theoretical  model  should  be  revised  in  this  case. 

Case  2:  Adiabatic  wail 

Figures  8a  and  8b  show  the  eigenvalue  distribution  of  the  second  and  first  modes  for  the 
adiabatic  plate,  T*  =  T'ad,  at  the  Mach  number  Me  =  6.8.  The  freestream  temperature  is 
T’  =  702578  K,  specific  heat  ratio  is  y  =  1.4 ,  Prandtl  number  is  Pr=0.72,  second 

viscosity  coefficient  is  k  =  1.2,  and  disturbance  frequency  parameter  is  F  =  1  O'4 .  It  is 
seen  that  the  functions  ar(R)  and  at(R)  are  smoother  compared  with  the  cooled  wall 
case  shown  in  Figs.  3a  and  3b.  This  is  due  to  the  fact  that  the  spectrum  branch  points  are 
shifted  from  the  real  axis  to  the  complex  R  -plane  as  shown  in  Fig.  11. 

Figures  9a  and  9b  show  the  second-mode  increment,  -a, .  =  -Im«2(i?) ,  and  the 
receptivity  functions  Gvib(R),  GS(R )  and  Ge(R)  respectively.  Because  the  first  branch 
point  Rh]  is  far  from  the  vibrator  center  point  R  (see  Point  2  in  Fig.  11),  receptivity 
peaks  are  smoothed  out.  Nevertheless  the  receptivity  function  maximums  are  still 
observed  near  the  lower  neutral  branch. 

Figures  10a,  10b,  10c  and  lOd  show  distributions  of  -a,(F),  Gvih(R),Gs(R)  and 
Ge  (R)  at  various  frequency  parameters  F.  Imaginary  part  of  the  branch  point  is  negative, 
Im(i?M)<0,  for  all  frequencies  considered  (see  Fig.  11).  As  the  Reynolds  number 
increases  along  the  real  axis,  the  first  branch  point  is  bypassed  from  the  upper  side.  The 
second  mode  (with  the  phase  speed  c  1  - 1  /  Me  as  R  decreases)  becomes  unstable. 
Because  the  second-mode  phase  speed  is  less  than  1,  there  is  no  discontinuity  of  the 
receptivity  functions  as  contrasted  to  the  cooled-wall  case  shown  in  Figures  5b,  5c  and  5d 
(Curves  4  and  5). 
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Similar  to  the  cooled  plate  case,  the  hypersonic  boundary  layer  is  most  sensitive  to 
vertical  velocity  perturbations,  if  they  are  in  resonance  with  boundary-layer  modes. 
Receptivity  to  the  surface  temperature  perturbation  is  relatively  low. 


Estimates  of  receptivity  effectiveness 

To  estimate  the  receptivity  effectiveness  we  consider  a  two-dimensional  vibrator  of  the 
shape 


J  0atx'<-//2,x'>//2 

{  a&t-l  I2<x'  <1 12 


(3.11) 


where  a  and  /  is  nondimensional  amplitude  and  streamwise  length  of  the  vibrator 
respectively;  their  dimensional  values  are  expressed  as  a  -  a-Jx*0  ve  /  U*e  and 

r  =  i  jxyju: :  is  dimensional  coordinate  of  the  vibrator  center,  which 

corresponds  to  x'  =  0 .  Fourier  transform  of  the  shape  function  (3.1 1)  is  expressed  in  the 
form 


+ro  ^ 

p(a)  =  fg(x')exp(-/ax')<ix'  = — sin(«//2). 
J  a 


(3.12) 


For  the  resonant  wavenumber  a  =  a0 ,  the  module  \p{a0)\  attains  its  maximum  value 


l/X«o)l=—  at  l  =  —(n  +  ^)  ,»  =  0,1,2, 

an  a n  V  2 J 


(3.13) 


Then,  in  accordance  with  Eq.  (3.4),  the  second-mode  wave  amplitude  is  evaluated  as 


«W,)  =  —  G,,„. 


(3.14) 
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For  the  adiabatic  flat  plate,  maximum  value  of  the  receptivity  function  is  Gvih  *  5  •  10  3  at 
the  frequency  parameter  F  =  10~4  and  Mach  number  Me  =  6.8  (see  Fig.  9b);  the 
resonance  wavenumber  is  a0  =  0.097  «  0.1  (see  Fig.  9a).  In  this  case,  the  vibrator  of 

amplitude  a *  »  O.OhJx'0ve  !U\  induces  the  second-mode  wave  of  initial  amplitude 
Qm(x0,a0)  ~  0.1% .  For  a  typical  hypersonic  wind  tunnel  with  the  unit  Reynolds  number 
Re,  =  107  1/m,  such  a  vibrator  has  the  amplitude  a’  =  10"6  m  and  length  /*  »  3  mm,  if 
it  is  located  near  the  lower  neutral  branch  point  x*  «  0.1  m. 

Note  that  the  highest  receptivity  occurs  near  the  low  neutral  branch.  Normal  waves 
generated  in  this  region  become  unstable  immediately  behind  the  force  element  and  can 
achieve  critical  amplitudes  earlier  than  those  generated  by  other  external  disturbances  far 
from  the  neutral  branch. 


4.  Receptivity  near  the  spectrum  branch  point 

In  Sections  2  and  3,  we  found  that  the  receptivity  function  tends  to  infinity  as  the 
resonance  point  tends  to  the  branch  point.  This  singularity  indicates  that  the  receptivity 
model  developed  in  Section  2  should  be  revised  in  the  branch  point  vicinity. 

In  this  Section,  we  analyze  the  distributed  excitation  of  the  second-mode  wave  by  force 
waves  when  the  resonance  point  tends  to  the  branch  point.  Then  we  consider  the 
instability  excitation  produced  by  a  local  force  element  placed  near  the  branch  point. 


4.1  Spectrum  behavior 

We  briefly  discuss  the  discrete  spectrum  behavior  in  the  branch-point  vicinity.  We 
assume  that  two  modes,  say  Mode  1  and  Mode  2,  coincide  at  the  branch  point  x,  =  xlb ; 
i.e.  al(xlb)  =  a2(xlb)  =  ab.  As  shown  in  [13]  and  [21],  such  branching  is  typical  for 
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hypersonic  boundary  layers,  when  the  first  mode  (Tollmien-Schlichting  wave)  and  the 
Mack  second  mode  coexist.  Near  the  branch-point,  eigenvalues  of  Modes  1  and  2  are 
expanded  as 


al2  =ab±ttjxl  -xlb+...,  (4.1a) 

Hereafter  the  subscript  b  denotes  quantities  at  the  first  branch  point  located  near  the 
lower  neutral  branch. 

Analysis  of  numerical  examples  shown  in  Figs.  3a,  3.b  and  8a,  8b  indicates  that 
|lm(A)|«|Re(A)|;  i.e.  Modes  1  and  2  are  almost  neutral  upstream  from  the  branch 

point.  For  simplicity  we  consider  the  case  when  the  parameters  A  and  xlb  are  real  and 
positive.  For  x,  <  x]b ,  the  eigenvalues  are  determined  as 

a[2  =ab  +^xlb  -x,  +...  (4.1b) 

Equation  (4.1b)  shows  that  Mode  1  has  smaller  wavenumber  (larger  phase  speed)  than 
Mode  2  in  the  upstream  region  x,  <  x16 .  If  the  branch  point  is  bypassed  from  the 
upper  side,  then 


al2  =ab±  i/L*Jx{  -  xlfe  + . . .  for  x,  >  xlfc  (4.1c) 

In  this  case,  Mode  1  is  stable  (lm(a,)>0)  and  Mode  2  is  unstable  (lm(«2)<0) 
downstream  from  the  branch  point. 

If  the  branch  point  is  bypassed  from  the  lower  side,  then 

al  2  =ab+  iA^jx]  -x]b  + . . .  for  x,  >  xlb  (4. 1  d) 
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In  this  case,  Mode  1  is  unstable  (lm(a2)<0)  and  Mode  2  is  stable  (lm(ai)>0) 
downstream  from  the  branch  point.  Note  that  in  both  cases  the  unstable  mode  is 
Mack’s  second  mode.  Without  loss  of  generality  we  consider  the  second  case 
assuming  that  Mode  1  is  unstable  downstream  from  the  branch  point  xxb  and  the 
eigenvalue  behavior  is  determined  by  Eqs.  (4.1a),  (4.1b)  and  (4. Id). 

As  the  streamwise  coordinate  x,  tends  to  the  branch  point,  Mode  1  strongly  interacts 
with  Mode  2  due  to  nonparallel  effects.  This  interaction  can  be  described  by  the  two¬ 
mode  approximation;  i.e.  the  disturbance  amplitude  from  (2.4)  is  expressed  in  the 
form 


F(xI,y)  =  c1(x1)4(xi,y)exp(ff'1/51)  +  c2(x1)22(x1,y)exp(f-,f52),  (4.2) 

<S|  2  =  2  (-X|  )fi?X|  . 

•V1  b 


As  shown  in  Ref.  [21],  the  amplitude  coefficients  c12(x,)  are  solutions  of  the  ODE 
problem 


=P-  =  cxWxx  +  c2W]2  exp [e~'i(S2  -  5, )] ,  (4.3a) 

dxx 

Hr 

-^--c2W22  +  cxW2X  exp [ff-'zXS,  -S2)],  (4.3b) 

dxx 

cMXs)  =  cXs,  c2(xXs)  =  c2s,  (4.3c) 

where  cXs  and  c2s  are  prescribed  values  at  a  certain  initial  point  xXs;  the  matrix 
elements  are  defined  as 


Wjk  = 


y  p  da  dx . 


i  +  i\  B j,Hx 


dx , 


B  8H 
1  da  dx. 


■,  y =1,2,  /t =1,2. 


(4.4) 
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Fedorov  and  Khokhlov  [13],  [21]  showed  that  the  matrix  elements  behave  as 


Wjk  x,  -> xlb . 


4(Xj  -xlfc) 


That  leads  to  the  following  singularity  of  the  amplitude  coefficients 


c,2(x,)  = 


C 


1,2 


(x,  -x,i) 


■  + . . .  as  x,  — >  x,h , 


(4.5) 


where  constants  C,  2  depend  on  the  normalization  of  eigenfunctions  Al2(x],y) .  Note 

that  this  singularity  is  totally  determined  by  the  spectrum  topology  prescribed  by  Eq. 
(4.1a).  In  the  first-order  approximation,  the  functions  cu(x, )  have  a  universal  form 

in  the  branch-point  vicinity;  they  do  not  depend  on  mean-flow  characteristics  and 
disturbance  parameters. 

In  the  region,  x,  -xib  -  0(s213 ) ,  the  mode  decomposition  (4.2)  is  not  valid. 
Introducing  the  inner  variable  £  =  e~2li(xl  -xlfc),  we  can  express  the  local 
disturbance  field  in  the  form 


F  =  s~'l6[c0(^)A0(y)  +  £~'l5c0^)A0X(y)  +  ...]ex^(i£-lSh  +i£~xnab<Z),  (4.6) 
where  c0(4")  is  solution  of  Airy  equation 

d2Cn  1 

—f~Jl2Cco=0.  (4.7) 


Fedorov  and  Khokhlov  [13],  [21]  matched  solutions  of  Eq.  (4.7)  with  the  normal¬ 
mode  expansions  of  the  outer  solution  and  established  an  exchange  rule  for  Modes  1 
and  2  in  the  branch-point  vicinity. 
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In  the  case  of  real  and  positive  X  and  xlb ,  a  structure  of  the  branch-point  influence 
domain  is  shown  in  Fig.  12.  It  includes  three  Stokes  lines  along  which  the  modes  are 
neutral.  Within  the  sector  containing  the  axis  portion  x{  >  xlb ,  Mode  1  is  unstable  and 
Mode  2  is  stable,  if  the  branch  point  is  bypassed  from  the  lower  side.  In  accordance 
with  (4.5)  and  (4.6),  the  disturbance  amplitude  has  a  peak  in  the  influence  domain, 
with  maximum  amplitude  being  of  the  order  of  s~U6 .  In  order  to  separate  this  effect 
from  the  receptivity  mechanism,  we  need  to  establish  a  relationship  between  the 
constants  C12  (instead  of  the  amplitude  coefficients  c12)  and  the  external  force 

characteristics. 


4.2  Distributed  receptivity 
Outer  solution 

We  consider  wall  perturbations  in  the  traveling-wave  form  given  by  Eq.  (2.2).  If  the 
resonance  point  x10  (where  the  force  wavenumber  a  coincides  with  the  unstable 
mode  eigenvalue  a0  =  Gr,(x10))  is  far  from  the  branch  point  xih,  then  the  normal 
wave  amplitude  F{  is  given  by  Eq.  (2.25).  In  the  absence  of  resonance  detuning 
(  /  =  0 ),  this  equation  is  expressed  as 


=  s~U2qd^Al(xi,y)ex p 


2m 


is 


Ja,(x,  )dX\ 


,  X]  >  X]Q  , 


(4.8) 


where  b  and  q  are  defined  by  Eqs.  (2.21b)  and  (2.21c)  respectively.  Differentiating 
Eq.  (2.10a)  with  respect  to  x,  and  imposing  the  solvability  condition  on  the  resulting 
equation  we  obtain 


u  da , 

b  =  -r±(xi0)  =  - 

flX, 


■®0  ’  -  (^10  »^o  )-^o 


dx. 


-  dH0  -  ' 

So.  a  (*|0>«oK 
I  da  , 


(4.9) 
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Using  this  relation  we  write  Eq.  (4.8)  as 


F,  =  s  'n*j2mb- 


(-/)(Q„,£0)V=0 

U  ZJh(x  „  w 

ox. 


-4(x,,.y)exp 


is 


Ja,(X[  )dX\ 


,  X,  >  X,0  . 


(4.10) 


If  the  resonance  point  tends  to  the  branch  point  as  x10  -*  x16  +  0,  the  eigenvalue 
derivative  b  is  expressed  from  Eq.  (4.1a)  in  the  form 


^  =  *t(*io  -xia)'1/2  +•• 


(4.11) 


Substituting  (4.1 1)  into  (4.10)  we  get  in  the  first  order  approximation 


i  » 


(4.12) 

Comparing  (4.12)  with  (4.2)  and  (4.5)  we  express  the  amplitude  coefficient  of  the 
unstable  Mode  1  in  the  form 


F  =s 


-1/2 


(®0’F0  )a=0 


■A{(xx,y)tx  p 


is~{  Jar,(x,; 


X[  >  x10 . 


Cx=s~xt24n^-j - ^°’^o)v=0 - r  as  xl0  ->x1A.  (4.13) 

{Bo,-j~(xlb,ab)A0\ 

Equation  (4.13)  shows  that  distributed  receptivity  has  no  singularity  at  the  spectrum 
branch  point.  In  accordance  with  Eq.  (4.5),  the  disturbance  peak  is  totally  due  to  a 
singular  behavior  of  the  unstable  mode  itself.  If  one  measures  the  disturbance 
amplitude  Q  at  a  certain  fixed  point  downstream  from  the  branch  point  for  various 
locations  of  the  resonant  point  x)0 ,  then  he  will  not  observe  a  peak  in  the  distribution 
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Q(xi0 )  at  the  point  x10  =  xlb .  However,  the  function  Q(xi0 )  will  have  a  maximum  at 
the  lower  neutral  branch  because  of  a  different  amplification  history  associated  with 
the  exponential  term  in  Eq.  (4.2).  Since  the  branch  point  xih  is  close  to  the  lower 

neutral  branch,  this  maximum  will  be  observed  at  xl0  «  xlh . 

Inner  solution 

We  assume  that  the  forcing  wavenumber  a  satisfies  the  resonance  condition  near  the 
branch  point;  i.e.  a  =  ab+£V3y.  Here  the  parameter  y  =  0(1)  characterizes 
resonance  detuning.  In  the  branch-point  vicinity,  the  disturbance  field  is  represented  in 
the  form 

F  =  [, s-2l3cb(C)Ab(y )  +  e'Uic0] (£)A0](y)  + ...  +  Ay(<Z,y)]exp(i£~U3abC) .  (4.14) 

Eigenfunction  Ah(y)  is  a  solution  of  the  homogeneous  problem  (2.10).  In  turn,  the 
function  A0]  (y)  is  a  solution  of  the  equation 

f  a  A  AfJ 

Am-—±(xlt,at)At=  0.  (4.15) 

yay  y  oa 

Substituting  the  expansion  (4.14)  into  Eqs.  (2.5)-(2.7)  we  obtain 

/  a  A  Ar  fill 

^-H0(xlb,ab)  c0lA0l+i^^(xlb,ab)Ab=  0.  (4.16) 

\dy  )  da 

From  Eqs.  (4.15)  and  (4.16)  we  derive  the  relationship 

e0,(O  =  -^-  (4-17) 

Expanding  the  matrix  H0  in  the  vicinity  of  xlb  we  obtain  the  following  problem  for 
the  forcing  term  Ay 
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(X  „  1  d2c„  r  1  d2H,  - 


0  -  rm9i 

(XIA,a6  )^40l  C*^T  O  (^li  )ab  )A  ~  ®  » 
m  <5c, 


( A )  1, 3.5.7  =  <P\x 3.4  ( )  exP(  *X )  at  j  =  0 , 


(4.18a) 


(4.18b) 


Av  — »  0  at  y  — »  cc . 


(4.18c) 


Using  a  standard  approach  (see  for  example  [28])  we  can  express  the  problem  (4.18) 
in  the  form 


~~H0(xlb,ab))Ar  +6bA(y)eiri  \^-^-{xxh,ab)Ah  + 

oy  )  dc,  2  <?er 


^-(xlb,ab)A0l  -cAxlb,ab)Ab--0, 
da  ax. 


(■^v )  1,3,5, 7  —  0  at 37  —  0, 


(4.19a) 


(4.19b) 


L5J  — >  0  at  37  ->  00 . 


(4.19c) 


Solution  of  this  problem  is  not  trivial,  if  the  inhomogeneous  term  of  Eq.  (4.19a)  is 
orthogonal  to  the  eigenfunction  Bb(y)  of  the  conjugate  problem.  This  solvability 
condition  leads  to  the  following  equation  for  the  amplitude  coefficient  ch{£) 


-L-A2£cb+qe'«  =  0, 


(4.20a) 


Ml 

1?H.*M1 


(4.20b) 


(Bb’^b  )  v=0 


1  U  d2Hn  - 


-  Bb,—^Ab  +  Bh,—±4 


2P’  7.  P*  ^  01 


r.  r?Hn 


(4.20c) 
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If  normal  waves  have  zero  amplitude  in  the  upstream  region,  i.e.  cb(£)-> 0  at 
£  — >  -oo,  then  solution  of  Eq.  (4.20a)  is  expressed  as 


cb  =-J%T  -  Bi(s)Ai(s')]exp(iyA~2l3s')ds' ,  s  =  A2I3£ ,  (4.21) 

^  -=0 

where  Ai(s)  and  Bi(s)  are  Airy  functions  [29].  For  the  branch-point  influence 
domain  structure  shown  in  Fig.  12,  we  obtain  the  following  asymptotic  expressions  in 
the  limit  s  -»  +co 


Ai(s)-> 


r~  „i/4 
y/KS 


exp 


f-ls 3/2 

V  3 


,  Bi(s)  -> 


r~  1/4 

yl  KS 


exp 


-s3/2 


4nq  (2  A 


^  ^-^rrl7rexPI 


,3/2 


7,  /=  J^4/(j)exp(zA  2nys)ds. 


(4.22a) 

(4.22b) 


Using  the  integral  form  of  Airy  function 


J  T-WJ 

Ai(s)  = —  )exp(/wi'  +  zw3  /3)du, 
2k  j 


(4.23) 


we  can  evaluate  (4.22b)  as 


7  = 


2k 


■*-00 

J<7z/exp(zw3  /3)  Jexp(zus  +  iX~2l2ys)ds  =  exp(-zX2;r3  /3).  (4.24) 


Summarizing  we  conclude  that  in  the  branch-point  vicinity,  x]  -x]b  =0(e213),  the 
force  wave  excites  the  unstable  normal  wave  of  the  amplitude 


Ft=s 


-1/2 


q^K 


Xi,2(x]  -x]b) 


-exp 


f  _  r3 
v"Z3 7*j 


4(*i,F)exp 


J/a,(x,  )dxt 


,(4.25) 


£  «X,  -  Xih  «1  . 
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From  Eqs.  (4.20b),  (4.20c)  we  get 


q  i-^b  >  )  v=o 


(4.26) 


Substituting  (4.26)  into  (4.25)  we  obtain  the  following  expression  for  the  amplitude 
coefficient  at  y  =  0 


Comparing  (4.13)  with  (4.27)  we  conclude  that,  in  the  first  order  approximation,  the 
inner  and  outer  solutions  lead  to  identical  expressions  for  the  amplitude  coefficient  of 
unstable  normal  wave. 

4.3  Local  receptivity 

Now  we  consider  a  local  force  element  with,  the  shape  function  determined  by  Eq. 
(2.26).  We  will  analyze  receptivity  as  a  function  of  the  distance  between  the  forcing 
element  and  the  branch  point  xlb . 

Outer  solution 

If  a  three-dimensional  force  element  is  far  from  the  branch  point  xib,  then  the 

amplitude  F]  of  unstable  wave-packet  generated  by  this  element  is  given  by  Eq. 

(2.32) .  Two-dimensional  forcing  element  excites  a  normal  wave  of  the  amplitude 

(2.33) .  Using  the  expression  (2.21c)  for  q,  the  relation  (4.9)  for  the  eigenvalue 
derivative  and  the  eigenvalue  expansion  (4.1a)  we  obtain  the  following  approximation 
for  q 


(4.27) 
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(®,AL 


Ap(a0) 


U  dH 0  -  \2(x10 -*,*) 

-®0>  a  (*10.««H  ) 

\  3*i  / 


1/2  aS  *•<>  ~ *  ■*!*  ' 


(4-28) 


Then  the  amplitude  coefficient  from  Eq.  (4.5)  is  expressed  in  the  branch  point  vicinity 
as 


C,= 


(O0,50)v=o 


Ap(a0) 


's  - \2(*10-x1A) 

*0,  (^'lO’^O  Mo  ) 

\  dxl  / 


1/4  as  x10  -»  xlfc 


(4.29) 


Local  receptivity  increases  proportionally  to  (xl0  -  xlb  )_l/4 ,  as  the  force  locus  tends  to 
the  branch  point.  This  effect  is  due  to  the  increase  of  the  force  wavenumber  range 
involved  into  the  excitation  process.  In  fact,  asymptotic  expansions  are  controlled  by 
the  parameter  sib.  If  the  forcing  element  is  far  from  the  branch  point,  then  the 
eigenvalue  derivative  b  =  0(1),  and  the  wavenumber  range  is  estimated  as 
{a  -  a0 )  =  0(sU2 ) .  As  the  forcing  element  tends  to  the  branch  point,  the  derivative  b 
tends  to  infinity  in  accordance  with  Eq.  (4.11).  If  the  force  element  is  placed  in  the 
region  ( xi0  -  x]b )  =  0(s213 ) ,  then  the  wavenumber  range  becomes  of  the  order  of 

(a  -  a0 )  =  0(sm ).  In  this  case,  Eq.  (4.29)  indicates  that  the  amplitude  coefficient  is 
of  the  order  of  C,  =  0(s~U6 ) . 

Inner  solution 

We  consider  a  force  element  placed  in  the  branch-point  influence  domain,  i.e. 
(*10  ~x\b  )  =  0(s2n ) .  According  to  Eq.  (4.25)  the  normal  wave  amplitude  generated 

by  the  forcing  wave  components  ab  +  sin(y  +  dy)  and  /?  +  d{3  is  expressed  as 


dF]  = 


-1/6 


(2  nY 


r\x{-xxb) 


-exp 


<  if} 

V  3A2 , 


exp 


iahs  +  —  ■ 
b  3 


2  A(x,  -x.b) 


p(ab  +  s'ny,j3)Al(xlb,y;/3)x 
eipz'dydp.  (4.30) 
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Integrating  over  y  and  (3  we  obtain  in  the  first  order  approximation 


F\  =£~m  |  ~~ --  --JA | (x, ,j;;/?)exp(z/fe')exp 
-Anyntf  (x,  -xlfe) 


where  the  integral  J  is  evaluated  as 


^ .  2  (x,  -x1A)3/2  V 

iabs  +  -  -  dp 

V  3  s  J 
(4.31) 


J=  \p(ab+£Uiy,P)ex p  -i^-\y  =  p{ab,P){bX2)n  Y{\n>) ,  (4.32) 

V  3/1  y  V3 


where  Gamma-function  is  f(l/ 3)  =2.678938.. 


In  the  case  of  a  two-dimensional  forcing,  the  integral  (4.31)  can  be  expressed  in  the 
compact  form 


_1/6  Xp{ab)T(m) 


2 4n{x{  -x1A) 


1/4 


X 


i^b^b  )v=0  ~A  .  . 

4(x,,j>)exp 


'sA 


*'l 

i  Jct^x,  )^x, 


V  *10 


<3*;, 


(4.33) 


where  £2/3  «x,  -x16  «1 .  From  this  expression  and  Eq.  (4.5)  we  obtain  the 
following  amplitude  coefficient 


C,= 


r£ V'* ipK)r(i/3)  (Bt,ot),„0 

1  3s  J  2V/r 


B  A 

Db’  -  A b 

\  OX, 


(4.34) 


Summarizing  we  conclude  that  local  receptivity  strongly  depends  on  a  distance 
between  the  forcing  element  and  the  branch  point.  As  xI0  ->  xib,  the  amplitude 

coefficient  C,  increases  proportionally  to  |x10  -xJ~l/4  (in  accordance  with  Eq. 
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(4.29))  and  attains  its  maximum  value  of  the  order  of  £~'16  at  the  point  xl0  =  xlh  (in 
accordance  with  Eq.  (4.34)).  If  one  measures  the  disturbance  amplitude  Q  at  a  certain 
fixed  point  x„  downstream  from  the  branch  point  for  various  locations  of  the  forcing 
element,  then  he  will  observe  a  local  peak  in  the  distribution  Q(xl0)  at  the  point 
*10  =*i b,  which  is  close  to  the  lower  neutral  branch.  This  peak  will  be  enhanced  by 
the  maximum  of  Q(xl0 )  due  to  the  amplification  history  in  the  region  x10  <  x,  <  x,„ 

associated  with  the  exponential  term  in  Eq.  (4.2).  Such  an  experiment  may  be 
conducted  on  a  cooled  flat  plate  at  the  Mach  number,  wall  temperature  ratio  and 
disturbance  frequency  corresponding  to  the  case  when  the  spectrum  branch  point  is 
close  to  the  real  axis.  For  example,  Fig.  4c  shows  such  a  regime  at  Me  =  5.95 , 

Tw  /  Tad  -  0. 1  and  frequency  parameter  F  =  1(T4  (see  Point  3). 

4.4  Results 

We  consider  excitation  of  Mack’s  second  mode  by  two-dimensional  local  forcing.  As 
discussed  in  Subsection  3.3,  this  mode  may  be  Mode  1  or  Mode  2  depending  on  the 
branch-point  bypass.  For  a  force  element  located  far  from  the  branch  point,  the 
receptivity  functions  Gvib ,  Gs,  and  Ge  are  determined  by  Eqs.  (3.5),  (3.7)  and  (3.9) 

respectively.  If  the  force  element  is  close  to  the  branch  point,  sin  «|x10  -x16[«l, 
then  the  amplitude  coefficient  is  determined  by  Eq.  (4.29). 

In  order  to  separate  the  receptivity  mechanism  from  the  mode  singularity  at  the  branch 
point,  we  normalize  the  receptivity  function  as 

_  ,  .j/4 

G  (xl0,xlb,ab )  =  G(x]0,xlb,(xb  )|x10  —  x16|  .  (4.35) 

Using  Eq.  (4.29)  we  obtain  the  following  expressions  for  the  receptivity  functions  in 
the  region  e213  « |x10  -  xlb  j « 1 : 
for  vibrations 
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^vib  (X]b’ab  )  ~ 


Qh(ym) 


(U  Bbl  +  icoBbi  +  T  Bb  5 )  v=0  x 


i  dx,  / 


for  suction-blowing  through  a  slot 


Gs(x]h,ab)  = 


Qb(ym) 


(^b.3  )  v=0 


X 


V^I 


for  surface  temperature  perturbation 


G*(*lfc,afc)  = 


(A.5  )  i  —O 


X 


,  (4.36a) 


(4.36b) 


(4.36c) 


Here  subscript  “6”  denotes  quantities  at  the  branch  point  where  Mode  1  and  Mode  2 
have  identical  eigenvalues  and  eigenfunctions. 

If  a  force  element  is  placed  at  the  branch  point,  xw  =  xl0 ,  then  we  get  from  Eq.  (4.34): 

for  vibrations 


Grib(xlh,ab)  = 


Qb(y,n’xib) 


(U'Bbi  +icoBbi+TBb5)y^Xih  x  f  ’  V'6 


'Mo 

i  dx , 


{xxb,ab)Ab,Bb 


r(i/3> 


A 

\$£j 


(4.37a) 


for  suction-blowing 


a  o.  ),_■  (^3>,~0 — r  -4=  ro  /3)f4. 


1/6 


,  (4.37b) 


for  surface  temperature  perturbation 


43 


f>ab )  ~ 


Qb(y,n)- 


(Bb, 5  )  v=0 


dE, 

dx. 


(x,b,ab)Ab,Bl 


7\ikrm(s 


,  1/6 


.  (4.37c) 


We  determine  the  global  length  scale  as  L  =  Re(x*b).  Then  the  branch-point 
coordinate  is  x16=l,  and  x,  =(R/Rb)2.  Receptivity  functions  GYjb,  Gs  and  Ge 
have  been  calculated  in  the  branch  point  vicinity  at  the  following  parameters:  local 
Mach  number  Me  =5.95 ,  flow  temperature  T*  =  70.2578  K,  wall  temperature  ratio 

rj  /  r;  =0.69539,  specific  heat  ratio  cp/cv=  1.4,  Prandtl  number  Pr=0.72,  and 

second  viscosity  coefficient  k  =  \2.  These  parameters  correspond  to  the  case 
considered  in  Section  3  and  illustrated  in  Figure  3. 

Figures  13a-13e  show  the  receptivity  function  Gv/i(x10 )  near  the  branch  point  xib  =  1 

for  the  frequency  parameter  F  x  104  =0.6,  0.8,  1.0,  1.2,  and  1.4  respectively.  Dotted 
line  shows  the  function  (?,.,■*  (x10 )  calculated  using  the  outer  solution  given  by  Eq. 
(3.5);  solid  line  with  square  symbols  shows  the  asymptotic  behavior  determined  by 
Eq.  (4.36a);  open  circle  shows  the  quantity  Grib(xl0  =xlh)  given  by  Eq.  (4.37a). 
Figures  13a  and  13b  correspond  to  the  case  when  the  branch  point  is  slightly  shifted  to 
the  lower  half  of  the  complex  x, -plain  (see  Points  1,  2  in  Fig.  6).  Figure  13c 
illustrates  the  case  when  the  branch  point  is  very  close  to  the  real  axis  (see  Point  3  in 
Fig.  6).  Figures  Id  and  le  show  the  case  when  the  branch  point  is  slightly  shifted  to 
the  upper  half  of  the  complex  x,  -plain  (see  Points  4,  5  in  Fig.  6). 

In  all  cases  under  consideration,  the  imaginary  part  of  the  branch  point  is  much 
smaller  than  its  real  part;  le.  the  real  axis  crosses  the  branch-point  influence  domain 
x,-x, b=0(e213)  schematically  shown  in  Fig.  12.  Because  the  parameter  2  is 

relatively  small  (its  value  is  about  2.5  xl0“2),  the  asymptotic  parameter  (2/f)l/6 , 
which  actually  controls  the  receptivity  maximum  given  by  Eq.  (4.37a),  is  not  very 
large.  It  is  ( A/£)u 6  «4.3  at  the  frequency  parameter  F'  =  6xl0"5  and  (A,/s)]l6  «3.2 
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at  F  - 1 .4  x  10  4 .  This  explains  relatively  low  value  of  Gvih  (xl0  =  xlh ) ,  which  restricts 

receptivity  peaks  at  the  branch  point.  Similar  reasoning  helps  to  explain  the  difference 
between  asymptotic  curves  (solid  lines  with  square  symbols)  and  actual  behavior  of 
the  outer  solution  (dotted  lines).  Nevertheless  our  asymptotic  model  captures  basic 
features  of  the  branch-point  effect  and  allows  us  to  estimate  receptivity  peaks  in  the 
branch-point  vicinity. 

Figure  14  shows  the  receptivity  maximum  G(xl0-xlb )  versus  the  frequency 
parameter  for  suction-blowing  through  a  two-dimensional  slot,  G  =Gs{xw  =  xlb); 
for  vibrations,  G  =  Gvib(xw  =  xlb ) ;  and  surface  temperature  perturbations, 
G  =  Ge  (x10  =  xlb ) .  Similar  to  the  case  of  regular  receptivity  shown  in  Fig.  9b  (where 

the  force  element  is  far  from  the  branch  point),  suction-blowing  is  the  most  effective 
source  of  instability  excitation  in  the  branch-point  vicinity.  Receptivity  to  temperature 
perturbations  is  essentially  lower  than  that  to  vibrations  and  suction-blowing.  Similar 
conclusions  can  be  made  using  distributions  of  the  receptivity  maximum  with  respect 
to  the  branch-point  locus  R e(Rb )  shown  in  Fig.  15. 


5.  Estimates  of  critical  vibrations  for  bypass 

Sufficiently  strong  forcing  is  able  to  generate  such  a  high  level  of  local  disturbances 
in  the  excitation  region  that  nonlinear  mechanisms  can  cause  early  breakdown  to 
turbulence.  In  this  case,  the  linear  stability  phase  is  bypassed  and  the  transition  onset 
occurs  just  behind  the  forcing  element.  To  estimate  the  force  level  providing  bypass 
we  assume  that  the  excitation  mechanism  is  linear  up  to  the  critical  amplitude  Qcr  of 
normal  waves. 

5.1  Local  vibrations 

We  consider  a  two-dimensional  vibrator  of  the  shape  given  by  Eq.  (3.1 1).  Its  Fourier 
transform  is  determined  from  Eq.  (3.12).  For  the  resonant  wavenumber  a  =  aQ,  the 
module  of  Fourier  component  p(a0)  is  given  by  Eq.  (3.13). 
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If  the  vibrator  is  far  from  the  branch  point,  then  we  can  use  the  results  of  Section  3.  In 
accordance  with  Eq.  (3.14),  the  vibrator  amplitude  providing  critical  level  of  mass- 
flux  disturbance  Qcr  in  the  excitation  region  (x,  -  x10 )  =  0{e 1/2 )  is  estimated  as 


o„(xl  0)« 


«0  Qcr 


(5.1) 


where  Gvib  is  determined  by  Eq.  (3.5). 

If  the  vibrator  is  placed  at  the  branch  point,  then  we  can  use  the  results  of  Section  4. 
The  critical  vibrator  amplitude  in  the  excitation  region  (x,  -  x]b)  =  0(s213 )  is 
estimated  as 


'  X 
kSJ 


abQcr 

2  Gvib 


where  Gvih  is  given  by  Eq.  (4.37a). 


(5.2) 


For  definiteness,  we  assume  that  bypass  occurs  at  Q  >  Qcr  =  1%.  Using  this  criterion 
we  calculated  the  critical  amplitude  acr  as  a  function  of  the  vibration  locus  x10 

expressed  in  terms  of  the  Reynolds  number  R  =  y]x*l0U*e  /  v*  . 


Figure  16  shows  the  distribution  acr(R)  for  adiabatic  flat  plate  at  Mach  number 
Me  =6.8;  mean-flow  parameters  correspond  to  the  case  shown  in  Figures  10a  and 

10b;  nondimensional  frequency  parameter  is  F  x 104  =  0.6,  0.8  and  1.4.  Since  the 
branch  point  is  relatively  far  from  the  real  axis,  calculations  were  conducted  using  Eq. 
(5.1).  Minimum  of  the  function  acr(R )  corresponds  to  the  region  of  maximum 
receptivity  (see  Fig.  10b).  It  is  close  to  the  lower  neutral  branch.  Minimum  value  of 
acr{R)  is  slowly  decreases  as  the  vibrator  frequency  decreases.  For  the  unit  Reynolds 

number  Re,  =  U*  lv*e  =106  1/ft  and  frequency  parameter  F  =  0.8  x  10“4,  the 
dimensional  critical  amplitude  is  estimated  as  acr  =  acr^v*ex*  !U]  «  0.001  inch. 
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Figure  17  shows  the  distribution  acr(Rh )  for  cooled  flat  plate  at  Mach  number 
Me  =5.95;  mean-flow  parameters  correspond  to  the  case  shown  in  Figures  14  and 
15.  Because  the  branch  point  is  close  to  the  real  axis,  calculations  were  conducted 
using  Eq.  (5.2).  Due  to  high  receptivity  in  the  branch-point  vicinity  the 
nondimensional  critical  amplitude  is  essentially  lower  than  in  the  previous  case.  For 
the  unit  Reynolds  number  Re,  =U*  lv*e  =106  1/ft  and  frequency  parameter 

F  =  0.8  x  10“4,  the  dimensional  critical  amplitude  is  estimated  as  a[r  « 0.0007  inch. 

These  examples  show  that  local  vibrations  of  the  resonance  frequency  are  effective 
sources  of  boundary  layer  tripping.  If  vibrations  are  equally  available  allover  the  plate 
surface,  initial  amplitudes  of  unstable  normal  waves  are  maximal  near  the  lower 
neutral  branch. 


5.2  Distributed  vibrations 

If  the  vibration  spectrum  contains  discrete  components,  then  the  receptivity  level  may 
be  essentially  higher  and  the  critical  amplitude  may  be  much  smaller  than  in  the  case 
of  local  vibrations.  This  situation  may  occur  when  one  of  discrete  modes  of  skin 
structure  vibrations  is  in  resonance  with  the  boundary-layer  modes. 

If  the  resonance  point  is  far  from  the  branch  point,  then  the  critical  amplitude  of  the 
vibrational  mode  is  evaluated  from  Eq.  (2.25)  as 


Qc 


rib  \ 


(± 

In 


1/2 


(5.3) 


where  b  and  Gvih  are  given  by  Eqs.  (2.21b)  and  (3.5)  respectively.  Using  Eq.  (4.9)  we 
can  express  the  critical  amplitude  in  the  explicit  form 


& 

Evib  \2nb 


.  1/2 


(5.4) 
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E\ib  (*K»«0  ) 


Qi(ym>x  io ) 


{U'B2X+i(oB2i+VB2'5)y^x^ 


■V10 


L  <&, 


(x,0  ,^Q  )-A2)B2 


If  the  resonance  point  coincides  with  the  branch  point,  then  the  mass-flux  disturbance 
amplitude  is  estimated  from  Eqs.  (4.14)  and  (4.21)  as 


Q~z-Via^Qi(ym,x\b)- 


(5.5) 


Using  Eq.  (4.26)  we  obtain 
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Then  the  critical  amplitude  of  the  vibrational  mode  is  estimated  as 


f-1 


2/3 


a 


(5.7) 


Figure  18  shows  the  distribution  acr{R)  evaluated  from  Eq.  (5.4)  for  adiabatic  flat 
plate  at  Mach  number  Me  =  6.8 ;  mean-flow  parameters  correspond  to  the  case  shown 

in  Figures  10a  and  10b;  the  frequency  parameter  is  T7  x  104  =  0.6.  0.8  and  1.4. 
Comparing  these  data  with  the  distributions  shown  in  Fig.  16  we  conclude  that  the 
critical  amplitude  of  the  vibrational  discrete  mode  is  much  less  than  that  of  the  local 
vibrator.  For  the  unit  Reynolds  number  Re,  =  U*e  tv]  =106  1/ft  and  frequency 
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parameter  F  =  0.8  x  1(T4,  the  dimensional  critical  amplitude  is  estimated  as  acr  = 


acrsjvlx *  IU*e  *  6.6x10  5  inch. 

Figure  19  shows  the  distribution  acr(Rb)  for  cooled  flat  plate  at  Mach  number 
Me  =5.95;  mean-flow  parameters  correspond  to  the  case  shown  in  Figures  14  and 
15.  Because  the  branch  point  is  close  to  the  real  axis,  calculations  were  conducted 
using  Eq.  (5.7).  Nondimensional  critical  amplitudes  are  essentially  lower  than  in  the 
previous  case.  For  the  unit  Reynolds  number  Re,  =  U*e  lv*e  =106  1/ft  and  frequency 
parameter  F  =  0.8xl0“4,  the  dimensional  critical  amplitude  is  acr  «9.3  x  10“6  inch. 

These  examples  indicate  that  discrete  vibrational  modes  of  the  resonance  frequency 
and  wavelength  are  extremely  effective  sources  for  the  boundary  layer  tripping. 


6.  Summary  discussion 

This  study  addresses  theoretical  modeling  of  hypersonic  boundary  layer  receptivity  to 
wall-induced  perturbations  on  a  flat  plate.  Since  the  receptivity  problem  depends  on 
different  scales  and  includes  small  parameters,  we  used  a  combination  of  asymptotic 
and  numerical  methods.  Asymptotic  analysis  allows  us  to  obtain  compact  analytical 
relations  for  the  receptivity  function,  identify  governing  nondimensional  parameters 
and  establish  basic  features  of  the  receptivity  mechanism.  Analytical  results  were 
integrated  in  a  computationally  non-intensive  module,  which  was  combined  with  a 
standard  stability  code.  This  module  provides  initial  amplitudes  of  normal  waves 
excited  by  wall  disturbances  of  various  types  and  shapes. 

If  the  discrete  spectrum  of  the  boundary-layer  disturbances  is  simple,  i.e.  eigenvalues 
of  normal  waves  are  different,  then  receptivity  to  wall  forcing  of  a  traveling  wave 
form  shows  the  following  features: 

•  Resonant  excitation  occurs  in  narrow  regions  of  the  length  l  =  Rb2S . 

•  Normal-wave  amplitude  is  proportional  to  {LIS)'11  =  R'/2 . 


49 


•  As  the  Reynolds  number  increases,  nonparallel  effects  (which  disturb  the 
resonance  condition)  decrease  and  the  relative  length  of  the  resonance  region, 
IIS,  increases.  In  the  parallel-flow  limit, RL  ->*>,  the  normal  wave  amplitude 
tends  to  infinity.  In  this  case,  the  upper  limit  of  disturbance  amplitude  is 
determined  by  nonlinear  effects. 

These  conclusions  are  consistent  with  the  results  of  [23]  obtained  for  subsonic 
boundary  layers. 

As  contrasted  to  the  case  of  distributed  forcing,  local  forcing  generates  the  boundary- 
layer  disturbance  of  a  finite  amplitude  even  in  the  parallel  flow  limit.  This  is  due  to 
the  fact  that  only  a  small  portion  of  the  forcing  spectrum,  (a-a0)  =  0(em),  is 
involved  into  the  resonance  mechanism.  As  s  ->  0 ,  the  resonant  excitation  increases 
proportionally  to  e~vl ,  whereas  the  spectrum  region  is  narrowed  down  proportionally 
to  em .  Because  both  trends  compensate  each  other,  the  normal-wave  amplitude  tends 
to  a  finite  value  as  s  -»  0 .  This  conclusion  is  consistent  with  the  results  of  [25] 
obtained  for  subsonic  boundary  layers. 

Calculations  showed  that  a  hypersonic  boundary  layer  is  very  sensitive  to  vertical 
velocity  perturbations  to  be  in  resonance  with  the  boundary-layer  normal  waves.  Such 
perturbations  may  be  induced  by  suction-bowing  through  a  hole  or  slot  as  well  as  by 
skin  vibrations.  Receptivity  to  surface  temperature  perturbations  is  relatively  low.  For 
adiabatic  flat  plate  at  Mach  number  Me  =  6.8 ,  a  two-dimensional  vibrator  of  the 

length  lxnla0,  amplitude  a  «  O.Ol^x'v’  !U\  and  frequency  F  =  10‘4  induces  the 

second-mode  wave  of  the  initial  mass-flux  amplitude  Qm(x0,a0)  «  0.1%  If  nonlinear 
breakdown  begins  at  the  critical  amplitude  Qcr  «  1%  ,  then  the  amplification  factor  is 
estimated  as  N  =  In (Qcr  /  Qin )  »  2.3 .  This  value  is  essentially  less  than  the  empirical 
value  TV  «  1 0  for  “quiet”  conditions. 

As  shown  in  [13]  and  [21],  discrete  spectrum  of  hypersonic  boundary  layer  has  branch 
points,  where  eigenvalues  of  the  first  and  second  modes  coincide.  In  the  branch-point 
influence  domain,  the  normal  wave  decomposition  is  not  valid  and  should  be  replaced 
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by  a  local  solution.  In  the  cooled  wall  case,  the  first  branch  point  is  located  very  close 
the  lower  neutral  branch  of  the  second  mode.  Both  propagation  and  excitation  of 
normal  waves  are  singular  in  the  branch-point  influence  domain.  The  receptivity 
function  tends  to  infinity  as  the  resonance  point  approaches  the  branch  point. 
Asymptotic  analysis  of  this  singularity  leads  to  the  following  conclusions: 

•  Distributed  receptivity  to  traveling  force  wave  is  not  singular  at  the  branch  point. 
The  receptivity  function  peak  is  totally  due  to  a  singular  behavior  of  normal 
waves.  If  one  measures  the  disturbance  amplitude  Q  at  a  fixed  point  downstream 

from  the  branch  point  xlb  for  various  loci  of  the  resonance  point  xl0 ,  then  he  will 
not  observe  a  peak  of  the  distribution  Q(x  10 )  at  the  point  xl0  =  xih . 

•  Local  receptivity  strongly  depends  on  a  distance  between  the  forcing  element 
locus  x10  and  the  branch  point.  As  xl0  -»  xlh ,  the  normal-wave  amplitude 

I  j-l/4 

increases  proportionally  to  |x10  -  x]b 1  and  attains  its  maximum  value  of  the 

order  of  e~U6  at  x10  =  xlh .  If  one  measures  the  disturbance  amplitude  Q  at  a  fixed 
point  downstream  from  the  branch  point  for  various  loci  of  the  forcing  element, 
then  he  will  observe  a  local  peak  in  the  distribution  Q(xl0 )  at  the  point  x10  =  xib . 

Suction-blowing  and  skin  vibrations  are  very  effective  sources  of  the  second  mode 
excitation  in  the  branch-point  vicinity.  Receptivity  to  surface  temperature 
perturbations  is  essentially  lower. 

Critical  amplitudes  acr  of  skin  vibrations  were  estimated  for  bypass  of  the  linear 
stability  phase.  Assuming  that  bypass  occurs  at  the  disturbance  level  Qcr  =1%,  we 
found  that 

•  for  a  local  two-dimensional  vibrator  on  adiabatic  wall  at  Mach  Me  =6.8,  the 
critical  amplitude  is  acr  «  0.09^/v‘x’  /U]  near  the  lower  neutral  branch; 

•  for  a  local  two-dimensional  vibrator  on  cooled  wall  at  Mach  Me  =  5.95 ,  the 
critical  amplitude  is  acr  «  0.025  fi*ex*  /  U]  near  the  first  branch  point; 
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•  for  two-dimensional  vibrational  wave  on  adiabatic  wall  at  Mach  Me-  6.8 , 
critical  amplitude  is  acr  * 0.005^/vJx*  IU*  near  the  lower  neutral  branch; 

•  for  two-dimensional  vibrational  wave  on  cooled  wall  at  Mach  M e  =  5.95 ,  critical 
amplitude  is  a„  «  0.0004^/v>*  IU'e  near  the  first  branch  point. 

In  all  cases  considered,  receptivity  maximum  corresponds  to  the  lower  neutral  branch. 
Normal  waves  generated  in  this  region  start  to  grow  just  behind  the  force  element  and 
can  achieve  critical  amplitudes  earlier  than  those  generated  elsewhere.  If  wall  forcing 
is  equally  available  allover  the  plate,  receptivity  near  the  lower  neutral  branch  is 

dominant. 

Foregoing  examples  show  that  hypersonic  boundaiy  layer  is  extremely  sensitive  to 
skin  vibrations.  This  factor  should  be  accounted  for  transition  prediction  in  active 
(ascending)  flight,  when  propulsion  system  may  induce  skin  vibrations  of  sufficiently 
high  levels.  This  conclusion  is  consistent  with  the  transition  data  [22]  obtained  on  a 
sharp  cone  installed  on  the  rocket  nose.  Transition  loci  measured  during  ascending 
flight  are  close  to  those  obtained  in  conventional  (noisy)  wind  tunnels. 

Our  theoretical  model  needs  to  be  verified  by  experiment.  This  can  be  done  in  a 
hypersonic  wind  tunnel  on  a  flat  plate  or  sharp  cone.  An  artificial  source  of  local 
vibrations  can  be  developed  using  high-frequency  piezoceramic  elements.  Another 
option  is  to  induce  periodic  suction-blowing  perturbations  through  a  hole  using  the 
spark-discharge  method  developed  in  the  Institute  of  Theoretical  and  Applied 
Mechanics  of  Novosibirsk  [30].  Our  theoretical  model  can  help  to  predict  optimal 
parameters  of  such  elements  and  estimate  disturbance  levels  to  be  measured. 

Of  special  interest  are  the  cooled  plate  regimes,  when  the  branch  point  effects  are 
most  pronounced.  These  effects  can  be  captured  by  hot-wire  measurements  of 
disturbance  field  generated  by  the  wall  forcing  element  of  a  fixed  frequency. 

Analytical  results  have  been  obtained  for  three-dimensional  disturbances.  However, 
the  numerical  examples  addressed  two-dimensional  cases  only.  In  the  majority  of 
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practical  cases,  disturbance  sources  are  three-dimensional.  Calculations  of  excitation 
and  evolution  of  unstable  wave  packets  including  the  branch-point  effects  may  give  a 
new  insight  into  the  initial  phases  of  hypersonic  boundary- layer  transition. 
Parametrical  studies  in  a  wide  range  of  Mach  number,  pressure  gradient,  wall 
temperature  ratio  and  other  parameters  are  also  needed. 
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Fig.  1  Flow  scheme. 
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Fig.  3b  Disturbance  spectrum:  Me  =  5.95,  Tw  /  Tad  =  0.1 ,  F  =  10~4 ;  distributions 
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Fig.  4d  Receptivity  function  Ge(R ) 
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Fig.  5a  Second-mode  increment  as  a  function  of  Reynolds  number:  Me  =  5.95 , 

Tw  /  Tad  =  0.1 ;  curves  1, 2,  3, 4  and  5  correspond  to  the  frequency  parameter 
Fx  104  =  1.4, 1.2,  1.0,  0.8  and  0.6  respectively. 
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Fig.  5b  Receptivity  function  Gvib(R)  at  Me  =  5.95 ,  Tw  /  Tad  =  0.1 ;  curves  1,  2,  3,  4 
and  5  correspond  to  F  x  104  =  1.4,  1.2,  1.0,  0.8  and  0.6  respectively. 
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Fig.  5d  Receptivity  function  Ge(R)  at  Me  =  5.95,  Tw  /  Tad  =  0.1 ;  curves  1,  2,  3, 4  and 
5  correspond  to  F  x  104  =  1.4, 1.2, 1.0,  0.8  and  0.6  respectively. 


Fig.  6  First  branch  point  locus  at  Me  =  5.95 ,  Tw  /  Tad  =  0.1 ;  points  1,  2,  3, 4  and  5 
correspond  to  F  x  104  =1.4,  1.2,  1.0,  0.8  and  0.6  respectively. 


Fig.  7a  Locus  of  the  receptivity  function  maximum  (solid  line)  and  the  lower  neutral 
branch  (symbols):  M  =  5.95 ,  T  /  T,  =  0.1 . 


Fig.  8a  Disturbance  spectrum:  Me  -  6.8 ,  Tw  /  Tad  =  1 ,  F  =  10  4 ;  distributions  ar(R ) 


Fig.  8b  Disturbance  spectrum:  Me  =  6.8 ,  Tw  /  Tad  =  1 ,  F  =  10  4 ;  distributions  a.  (F) 


Fig.  9b  Receptivity  functions  at  Me  =  6.8 
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Fig.  10c  Receptivity  function  GS(R )  at  Me  =6. 8,  Tw  /  Tad  =  1 ;  curves  1,  2,  3, 
and  4  correspond  to  F  x  104  =  1 .4,  1 .0,  0.8  and  0.6  respectively. 


Fig.  lOd  Receptivity  function  Ge(R)  at  Me  =  6.8 .  Tw/  Tad  -  curves  1,  2,  3, 
and  4  correspond  to  F  x  104  =1.4, 1.0,  0.8  and  0.6  respectively. 


Fig.  1 1  Branch  point  locus  at  Me  =  6.8  »  /  Tad  -  1  :  curves  1, 2,  3,  and  4  correspond 

to  the  frequency  parameter  F  x  104  =1.4,  1.0,  0.8  and  0.6  respectively. 


Fig.  12  Structure  of  the  branch-point  influence  domain. 


Fig.  13a  Receptivity  function  Gvib(x]0 )  near  the  branch  point  xlb  =  1 ,  F  =  0.6  x  1 0"4 ; 
dotted  line  -  outer  solution  (3.5),  solid  line  with  square  symbols  -  asymptotic 
approximation  (4.36a),  open  circle  -  Gvib(x]0  =xlb )  from  (4.37a). 


Fig.  13b  Receptivity  function  Gyib(xw)  near  the  branch  point  xlb  =1,  F  =  0.8xl0'4; 
dotted  line  -  outer  solution  (3.5),  solid  line  with  square  symbols  -  asymptotic 
approximation  (4.36a),  open  circle  -  Grib(xl0  =  xlb )  from  (4.37a). 
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Fig.  13c  Receptivity  function  Gvib  (x10 )  near  the  branch  point  x]h  = 1 ,  F  = 10“4 ; 

dotted  line  -  outer  solution  (3.5),  solid  line  with  square  symbols  -  asymptotic 
approximation  (4.36a),  open  circle  -  Grib(x  10  =  xlh )  from  (4.37a). 


Fig.  13d  Receptivity  function  Grib  (x10 )  near  the  branch  point  x]b  =  1 ,  F  =  1.2  x  10  4 ; 
dotted  line  -  outer  solution  (3.5),  solid  line  with  square  symbols  -  asymptotic 
approximation  (4.36a),  open  circle  -  Gvib(xXQ  =  xxb )  from  (4.37a). 
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Fig.  1 3e  Receptivity  function  Gvih  (x10 )  near  the  branch  point  xlb  =  1 ,  F  =  1 .4  x  1  O'4 ; 

dotted  line  -  outer  solution  (3.5),  solid  line  with  square  symbols  -  asymptotic 
approximation  (4.36a),  open  circle  -  Gvib(xi0  =  x]h )  from  (4.37a). 
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Fig.  14  Receptivity  maximum  G(xl0  =  xlh )  versus  frequency  parameter. 
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Fig.  15  Receptivity  maximum  G(xl0  =*„)  versus  the  branch-point  locus  Ret  ft,  ) 
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Fig.  16  Critical  amplitude  of  local  vibrations:  Me  =  6.8,  Tw  /  Tad  =  1 . 
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Fig.  17  Critical  amplitude  of  local  vibrations  near  the  branch  point: 
Me=  5.95,  TJTad=  0.1. 
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Fig.  18  Critical  amplitude  of  distributed  vibrations:  Me  =  6.8 ,  Tw  /  Tc 
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Fig.  19  Critical  amplitude  of  distributed  vibrations  near  the  branch  point:  Me  =  5.95, 


TJTad=  0.1. 


